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Abstract 

Preliminary  measurements  are  described  of  electron  number  density  obtained  in  the  near  field  of  a 
200W  Busek  Hall  thruster.  The  approach  taken  is  ultra-high  high  frequency  (3.3  mm  wavelength,  90 
GHz)  microwave  interferometry  that  affords  high  spatial  resolution  suitable  for  studying  the  near  exit 
region,  where  the  plasma  density  can  be  as  high  as  1 01  'em'3.  A  90  GHz  beam  can  be  focused  down  to  a 
waist  of  approximately  7  mm,  and  over  a  50  mm  plasma  path  length,  a  10ncm'3  plasma  density  gives 
rise  to  an  easily  measurable  3°  phase  shift.  The  system  is  of  a  phase-bridge  design,  utilizing  two  signal 
anus  split  from  a  fixed  frequency  source  (one  passing  through  the  plasma)  that  recombine  at  two 
balanced  mixers.  A  line-integrated  electron  density  is  obtained  by  comparing  the  two  signals  from  the 
mixers.  This  interferometer  is  suitable  for  measuring  time-dependent  plasma  density  fluctuations 
offering  unprecedented  information  about  plasma  turbulence  in  the  near  exit  region,  and  the  opportunity 
to  study  turbulence-enhanced  electron  transport  in  regions  where  the  plasma  is  presumably  collisionless 
and  free  of  interactions  with  the  channel  wall. 


I.  Introduction 

One  of  the  most  important  properties  to 
characterize  in  a  Hall  thruster  discharge  is  the 
near-field  plasma  (electron  or  ion)  density.* 
Together  with  the  ion  velocity,  the  plasma 
density  provides  an  independent  measurement  of 
the  near-field  ion  current.  The  local  plasma 
density  also  serves  as  a  benchmark  for  thruster 
simulations,  as  its  established  value  is  sensitive 
to  both  the  ionization  process,  and  the 
accelerating  electric  field.  The  most  commonly 
used  methods  of  determining  the  ion  density  in 
the  near-field  of  Hall  thruster  discharges  is 
through  the  use  of  either  cylindrical  Langmuir 

*  Hall  thrusters  are  quasi-neutral  devices,  so  the 
electron  and  ion  densities  are  to  a  great  extent  equal. 
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probes  in  current  saturation  mode  [1,2],  or 
guarded  ion  probes  [3],  both  of  which  suffer 
from  being  intrusive,  and  from  uncertainties  in 
their  interpretation  due  to  ion-enhanced 
secondary  electron  emission.  A  measurement  of 
time-average  plasma  density  using  probes  in  the 
near-field  of  Hall  thrusters  is  at  best  accurate  to 
an  order  of  magnitude,  and  time-resolved 
measurements  must  account  for  limited  spectral 
response  of  the  probe  construction.  On  the  other 
hand,  the  ease  of  interpretation  and  non-intrusive 
nature  of  alternative  diagnostics  such  as 
microwave  interferometry  makes  it  an  attractive 
choice  for  collecting  electron  (plasma)  density 
data.  Microwave  interferometry  at  centimeter- 
scale  wavelengths  (e.g.,  20  GHz)  has  been  used 
to  study  the  distant  plume  region  of  a  Hall 
thruster.  However,  its  use  in  studying  the  near¬ 
field  is  hampered  by  the  limited  spatial 
resolution  that  can  be  obtained  because  the  beam 
waist  when  focused  is  at  best  a  few  wavelengths 
in  diameter. 

In  this  paper,  we  present  preliminary  results 
of  measurements  of  path-integrated  plasma  line 
density  in  the  near-field  of  a  BHT-200  Hall 
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thruster  [4],  A  custom-designed  phase-bridge 
interferometer  operating  with  a  3 -mm 
wavelength  solid-state  microwave  source 
resulted  in  modest  spatial  resolution  (~7  mm) 
and  lateral  and  axial  surveys  of  the  plasma  line- 
density  within  50  mm  of  the  exit  plane  of  the 
thruster.  We  also  provide  a  brief  discussion  of 
some  of  the  challenges  encountered  with  drifts 
in  the  relative  sensitivity  and  phase  shift 
introduced  by  the  mixers,  due  to  heat  loading  of 
the  microwave  components  as  the  entire 
interferometer  is  operated  in  vacuum.  These 
drifts  precluded  the  measurement  of  plasma 
density  to  no  better  than  a  factor  of  about  two  to 
four,  and  were  the  main  source  of  experimental 
uncertainty.  Because  of  the  interferometer 
design,  measurements  can  be  obtained  of  the 
time-variation  in  the  plasma  density,  and  a  first 
direct  comparison  is  made  in  the  non-intrusive 
plasma  density  fluctuations,  to  fluctuations  in 
overall  discharge  current. 

II.  Background  Theory  and  Analysis 

For  collisionless  plasmas,  the  line-integrated 
electron  density  can  be  expressed  in  terms  of  the 
measured  phase  difference  between  two 
identical  microwave  sources,  one  of  which  is 
passed  through  the  plasma  of  interest  (signal 
arm),  the  other  serving  as  a  phase  reference 
(reference  arm). 

The  phase  difference  between  these  two 
waves,  introduced  by  the  complex  impedance  of 
the  free  electrons,  is  related  to  the  line-averaged 
electron  density  as  [5]: 


n 


e 


Anme  eQc2 
e2Xal 


A  <j) 
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Here  me  is  the  electron  mass,  e  is  the  electron 
charge,  c  is  the  speed  of  light,  X„  is  the  free  space 
microwave  wavelength,  e0  is  the  free-space 
permittivity  and  l  is  the  plasma  column  width. 
Although  an  interferometry  measurement  of  this 
sort  is  averaged  along  the  path  of  microwave 
propagation  though  the  plasma,  Hall  thrusters 
are  reasonably  axisymmetric,  so  that  radial 
variations  in  the  plasma  density  can  be 
reconstructed  from  a  series  of  measurements 
made  along  paths  passing  through  various 


chords  of  the  plume,  using  standard  inverse 
transform  methods  (e.g.,  Abel  transformations) 
or  numerical  “onion  peel”  techniques  such  as 
those  commonly  used  in  optical  absorption 
spectroscopy  for  characterizing  symmetric 
atomic  sources.  In  a  ground  test  facility,  the  Hall 
Thruster  can  be  translated  in  the  axial  and  radial 
dimensions,  allowing  the  microwave  beam  to 
pass  through  many  sections  of  the  plasma 
column  facilitating  this  inversion  process.  The 
inversion  process,  which  begins  from  the  outer¬ 
most  region  of  the  plasma  towards  the  center,  is 
most  accurate  where  lateral  variations  in  the 
phase-shift  are  the  greatest.  In  this  paper,  no 
attempt  has  yet  been  made  to  carry  out  this 
inversion  process. 

A  schematic  (block)  diagram  of  the  phase 
bridge-interferometer  configuration  used  in  this 
study  is  shown  in  Fig.  1  below.  In  the 
experimental  set-up,  the  output  from  a  90  GHz 
fixed  frequency  Gunn  diode  with  integral 
oscillator  (Millitek  model  GDM-10-0017IR) 
having  a  50mW  maximum  output  power  is  split 
into  two  arms  of  equal  strength  through  the  use 
of  a  high  directional  coupler  (Millitek  model 
CSS-10-R090BZ).  One  arm  leads  to  two 
variable  phase  shifters  (Millitek  models  VPS-10- 
R000N  providing  a  combined  360°  phase  shift) 
and  level  set  attenuator  (Millitek  model  LSA-10- 
R000N).  The  phase  shifters  and  attenuator  serve 
as  a  means  of  calibrating  the  interferometer.  The 
beam  is  then  launched  into  the  plasma  through  a 
standard  pyramidal  gain  horn  (Millitek  Model 
SGH-10-RP000)  and  custom  designed  matched 
Teflon  lens  (50  cm  focal  length)  with  an 
estimated  beam  waist  of  ~7  mm  diameter.  The 
beam  is  subsequently  collected  with  a  second 
matched  lens/hom  combination.  Both  the  return 
arm  from  the  collection  horn,  and  the  second 
arm  originating  at  the  diode,  are  split  further  into 
two  equal  sources  with  two  more  highly 
directional  couplers  (Millitek  model  CSS-1 0- 
R090BO).  The  two  pairs  of  source  (incident) 
and  return  arms  are  combined  in  two 
independently  balanced  and  biased  mixers 
(Millitek  model  MXB-10-RR0WX).  The  mixers 
serve  as  wave  multipliers,  with  low-pass  filters, 
having  a  DC-510  MHz  IF  frequency  response. 
All  components  of  the  microwave  circuit  are 
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coupled  through  the  use  of  sections  of  non¬ 
standard  length  WR-10  waveguides. 

For  an  incident  wave  S,(t)  =  Ajcos(cot)  and 
return  wave  Sft)  =  A2cos(cot+tj)),  the  mixer 
output  gives  a  DC  signal  proportional  to  the 
phase  shift  A <f>  introduced  by  the  plasma: 

/  =  cos  (f)  @) 

2 

The  source  from  one  arm  is  passed  through  a  90° 
phase  shifter.  For  the  branch  that  includes  90° 
phase  shifter,  the  mixer  signal  is: 

_  A,A2  (n  A,A2  .  /±\ 

Q  =  -A-A-cos  --(j)  =-~ sin ((f))  (3) 
2  \2  )  2 


Figure  1.  Schematic  diagram  of  the  phase  bridge- 
interferometer.  1-90  GHz  Gunn  Diode,  2- 
Splitters,  3-  Balanced  Mixers,  4  -  Variable 
Attenuator,  and  5  -  Variable  Phase  Shifter. 


We  see  then,  that  the  ratio  of  the  two  signals 
generated  in  the  balanced  mixers,  R  =  I/O  -  tan 
(Atpj  por  smau  phase  shifts,  where  tan  )  ~ 

A0 ,  then  this  ratio  is  ideally  equal  to  the  phase 


shift,  since  the  signal  amplitudes  cancel.  In 
principal,  fluctuations  in  the  output  of  any  single 
mixer  signal  is  proportional  to  the  fluctuations  in 
the  plasma  density.  The  ratio  of  the  mixer 
outputs  directly  reflects  fluctuations  in  the 
plasma  density. 

In  practice,  the  combined  branches  that  give 
rise  to  the  I  and  Q  signals  have  unequal  intrinsic 
losses  (for  example,  due  to  differing  waveguide 
properties  and  lengths,  amongst  others). 
Furthermore,  the  two  mixers  themselves  may 
have  a  different  sensitivity  to  signals  at  their 
inputs  (tti,  ocq)  and  an  intrinsic  DC  offset  (/„, 
Qo).  As  a  result,  the  mixer  signals  can  be 
represented  as: 


r  t  ^1^2  a  j. 

I  =  f0+aj  -  ■  --  cos A 0 


Q  -  Qo  +  ac 


^-sin  A  (j) 


Tire  signals  /  and  Q  trace  an  ellipse  in  Cartaesian 
space  (so-called  Smith  circles),  the  center  of 
which  is  located  at  position  (/„,  Q0),  as  illustrated 
in  Fig.  2.  In  general,  a  non-zero  reference  (no 
plasma)  phase  shift  <f>0  can  arise  due  to  intrinsic 
phase  shifts  between  the  two  paths.  In  Figure  2, 
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Figure  2.  Smith  diagram  of  the  microwave 
interferometer  response,  obtained  by  varying  the 
phase  setting  in  the  Phase  Shifters. 

the  calibration  (Smith)  circle  is  shown  for 
successive  5-degree  settings  in  the  variable 
phase  shifter.  Tn  this  figure,  the  calibration 
circles  have  not  been  re-scaled  for  phase  shift 
and  DC-offsets.  It  is  noteworthy  that  the  plot  is 
not  circular,  indicative  of  the  difference  in  the 
sensitivity  of  the  two  mixers.  Figure  3  provides 
a  series  of  Smith  circles  obtained  while  varying 
the  attenuators  by  5  dB  intervals.  In  this  figure, 
the  calibration  circles  have  been  re-scaled  for 
DC  offsets  and  intrinsic  phase  shift  (to  the  0  dB 
circle).  Note  that  although  «  ~1  dB  of 
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attenuation  is  expected  for  the  collisionless 
plume  of  this  discharge,  these  calibration  circles 
indicate  that  the  variable  attenuator  introduces 
an  added  phase  shift  which  would  have  to  be 
corrected  for  when  probing  collisional  plasmas. 
All  analysis  presented  here  use  zero-attenuation 
calibration  circles  for  data  reduction. 

It  is  noteworthy  that  with  the  interferometer 
located  within  the  vacuum  chamber,  and  the 
intrinsic  heat  dissipation  from  the  powered 
mixers  and  diodes  (as  well  as  some  heat  loading 
from  the  plume/thruster),  noticeable  drifts  were 
seen  in  the  sensitivity,  phase,  and  offsets  in  the 
mixers.  These  drifts  precluded  the 
measurements  of  plasma  density  to  no  better 
than  a  factor  of  two-five.  Since  making  the 
measurements  discussed  here,  we  have  made 
provisions  for  active  cooling  of  the  components, 
and  improved  accuracy  and  sensitivity  is 
expected  from  future  measurements. 


Q-Mixer  Signal  (V) 


Figure  3.  Smith  diagram  of  the  microwave 
interferometer  response,  obtained  by  varying  the 
phase  setting  in  the  Phase  Shifters  and  attenuation 
settings  in  the  Variable  Attenuator. 

Because  of  possible  drifts  in  the  zero-phase 
reference,  our  analysis  could  not  rely  on  the 
identification  of  the  precise  location  of  the  /  and 
Q  signals  on  the  calibrated  Smith  circles  when 
the  microwave  beam  passes  through  the  plasma 
plume.  Instead,  an  alternate  approach  to  data 
interpretation  is  taken  here.  The  analysis  of  the 
data  for  the  results  presented  here  is  relatively 
straightforward,  and  relies  on  identifying 
(experimentally),  a  zero-phase  shift  reference 


point  by  recording  mixer  signals  ( Iop ,  Qop)  with 
the  microwave  beam  passing  through  regions  of 
the  plume  where  there  is  little/no  plasma.  For  the 
small  phase  shifts  anticipated  when  passing 
through  the  near-field  plume  (<  5°),  the  arc 
swept  along  the  Smith  circle  from  this  zero- 
phase  shift  condition,  S  =  [(1-1  op)2  +  (Q-0„p)2]''2 
=  .  Here,  M  is  the  modulus  in  the  mixer 

signal  at  zero-phase,  corrected  for  the  DC  offset: 
AT  =  ((IoP-lo)-+(Qop-Qo)2)1/2.  The  phase  shift 
(and  hence  plasma  density  from  Equation  1)  is 
determined  then,  from: 


(/-U+b-g JT 


(4) 


A  representative  scan  of  the  l-I0Pi  and  Q-Q„P 
signals  obtained  as  the  thruster  is  traversed 
across  the  microwave  beam  at  an  axial  position 
of  z  =  95  mm  (z  =  0  mm  is  the  location  of  the 
exit  plane)  is  shown  in  Fig.  4. 
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Figure  4.  Mixer  signals  obtained  as  thruster  is 
traversed  across  microwave  beam. 

The  Resulting  lateral  variation  in  signal  S  is 
shown  in  Fig.  5.  An  estimate  of  the  modulus,  AT, 
for  the  mixer  signal  at  zero  phase  in  the  absence 
of  the  plasma,  obtained  from  a  calibration  circle 
carried  out  prior  to  the  measurements,  is  M~  60 
mV.  The  measurement  of  S  along  with  AT  is  used 
to  determine  the  phase  shift,  which  would  have  a 
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Figure  5.  Signal,  S  corresponding  to  the  arc  swept 
along  the  Smith  circle  as  the  plume  traverses  the 
microwave  beam. 

lateral  distribution  similar  in  shape  to  the  signal 
depicted  in  Fig.  5. 

III.  Experiment 

A  photograph  of  the  interferometer 
positioned  within  the  vacuum  is  shown  in  Fig.  6. 
To  aid  in  the  visualization  of  the  interferometer 
configuration,  a  three-dimensional  illustration  of 
the  interferometer  assembly  as  constructed  on  a 
rigid-frame  structure  is  shown  in  Fig.  7. 

The  Flail  thruster  studied  here  belongs  to  a 
family  of  Busek  BFIT-200  thrusters,  which 
normally  run  at  a  power  of  approximately 
200W.  It  has  an  inner  wall  diameter  of  14  mm, 
and  an  outer  wall  diameter  of  31  mm  (8  mm 
channel  width).  The  central  pole  piece  extended 
approximately  7  mm  beyond  the  exit  plane  of 
the  channel.  The  channel  exit  served  as  the 
origin  position  marker  for  the  axial  ( z )  position 
referred  to  in  the  measurements.  The  cathode 
was  located  approximately  20  mm  downstream 
of  the  central  pole  piece. 

Measurements  were  carried  out  in  a  chamber 
1 .8  m  in  diameter  and  3  m  long  with  a  measured 
pimping  speed  of  -32,000  1/s  on  xenon, 
provided  by  four  single-stage  cryo  panels  and 
one  50  cm  diameter  two-stage  cryo  pump.  The 
thruster  was  operated  at  nominal  conditions  of 
0.8  A  of  discharge  current,  250V  discharge 
voltage  (controlled  mode),  an  anode  mass  flow 


Figure  6.  Photograph  of  interferometer  installed  in 
the  vacuum  chamber. 


Figure  7.  Schematic  of  the  phase  -  bridge 
interferometer.  1-90  GHz  Gunn  Diode,  2-  Splitters, 
3-  Balanced  Mixers,  4  -  Variable  Attenuator,  and  5  - 
Variable  Phase  Shifter. 

rate  of  8.5  seem,  and  a  cathode  mass  flow  rate  of 
1  seem. 

IV.  Results 

Measurements  were  made  at  nominal 
discharge  operating  conditions  while  translating 
the  FTall  thruster  through  a  range  of  axial  (z)  and 
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lateral  positions  so  that  the  focused  microwave 
beam  passes  across  a  chord  of  the  symmetric 
plume,  in  a  plane  nonnal  to  the  z-direction. 
Mixer  signals  were  recorded  digitally,  and  the 
signal  modulus  S  and  reference  modulus  M  was 
then  used  to  determine  the  phase  shift  and 
plasma  line-density, «  l  .  The  variation  in  this 

line-density  with  lateral  position  is  shown  in  Fig. 
8. 

The  results  shown  in  Fig.  8  have  not  yet  been 
deconvolved  through  a  radial  inversion  process 
even  though  the  deconvolution  (Abel  inversion) 
process  is  greatly  simplified  if  we  can  assume 
that  the  plume  is  symmetric  about  the  Hall 
thruster  axis.  If  we  were  to  use  the  typical  width 
of  the  distribution  shown  as  a  path  length  (~  50 
mm),  then  the  average  plasma  density  across  this 
50  mm  path  length  is  about  1.6  x  1011  cm-3, 
consistent  with  what  is  expected  in  the  very  near 
field  (within  one  engine  diameter)  of  this  Hall 
thruster,  and  with  the  probe  measurements  of 
Beale,  et  al  [6],  We  have  found  that  in  the  very 
near  field  (e.g.,  45  mm),  interferences  with  the 
microwave  beam  are  introduced  by  the 
placement  of  the  cathode  (positive  lateral 
positions  in  the  figure),  and  so  data  was  obtained 
only  for  one  side  of  the  thruster.  It  is  noteworthy 
that  at  larger  axial  positions,  there  is  an  apparent 


asymmetry  which  we  attribute  to  the  cathode  jet. 

Measurements  were  also  taken  at  a  fixed 
lateral  position  corresponding  to  chords  passing 
through  the  axis  of  the  thruster  and  varying  axial 
positions.  The  variation  in  the  plasma  line- 
density  with  axial  position  is  shown  in  Fig.  9. 
The  line  density  is  seen  to  fall  exponentially 
with  axial  position  from  the  exit  plane  of  the 
thruster,  as  expected  for  a  diverging  plume.  The 
divergence  of  the  plume  is  a  consequence  of  the 
local  electrostatic  fields,  shaped  mainly  by  the 
applied  magnetic  field. 

This  phase-bridge  interferometer  also 
provides  a  non-intrusive  means  of  characterizing 
temporal  fluctuations  in  the  plasma  density.  For 
small  phase  shifts,  the  output  from  a  single 
mixer  will  have  a  temporal  component  that  is 
proportional  to  the  fluctuations  in  the  plasma 
line-density.  In  the  vicinity  of  the  exit  plane, 
large-scale  fluctuations  are  attributed  to  the  so- 
called  “breathing”  mode  of  Hall  discharges  [7], 
which  arises  as  a  result  of  the  competition 
between  rapid  depletion  of  mass  within  the 
channel  due  to  ion  acceleration,  and  the 
replenishing  of  xenon  through  anode  injection. 
This  large  scale  fluctuation  is  usually  identified 
through  analysis  of  the  temporal  fluctuation  in 
the  discharge  current.  However,  such  a  “global” 
parameter  may  not  necessarily  reflect  the 
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Figure  8.  Typical  results  of  plasma  line  density  versus  axial  position  obtained  from  the  interferometer.  The  exit 
plane  is  at  z  =  0  mm.  These  data  have  not  yet  been  inverted. 

6 


O  z  =  75  mm 
□  z  =  95  mm 
_  i  z  =  85  mm 


Approved  for  public  release;  distribution  unlimited. 


2004  AIAA  Joint  Propulsion  Conference,  Ft.  Lauderdale,  FL  July  11-14,  2004 


AIAA-2004-3775 


20  40  60  80  100  120  140 


Axial  Position  (mm) 

Figure  9.  Axial  variation  in  the  plasma  line-density  in 
the  Hall  discharge  plume. 

behavior  of  the  discharge  at  a  particular  position, 
such  as  the  exit  plane  or  near-field  of  the 
thruster. 

Figure  10  below  compares  the  Fourier 
spectrum  of  the  current  fluctuations  to  the 
fluctuation  in  plasma  density  seen  along  a 
central  chord  at  an  axial  position  of  z  =  45  mm 
from  the  exit  plane.  Apparent  in  both  traces  is 
the  relatively  strong  principal  oscillation  at 
about  17  kHz,  along  with  four  harmonics.  There 
is  a  remarkable  similarity  between  the  two  traces 
depicted  in  the  figure,  which  implies  that  the 
low-frequency  large  scale  oscillations  in  the 
plume  plasma  density  (near-field)  is  indeed 
representative  of  the  breathing  mode  that  is 
ubiquitous  in  the  low  frequency  oscillations  in 
discharge  current. 

The  microwave  diagnostic  permits  the  non- 
intrusive  recording  of  high  frequency  (»  100 
kHz)  disturbances  in  the  near  field  of  the 
thruster,  which,  prior  to  this,  have  been 
characterized  nearly  exclusively  by  intrusive 
probes.  The  only  drawback  associated  with  the 
application  of  microwave  interferometry  to 
study  these  fluctuations  is  that  it  is  a  line-of- 
sight  diagnostic,  and  will  resolve  only  those 
disturbances  that  are  of  relatively  large  scale  - 
possibly  azimuthal  waves  limited  to  mode 
numbers  less  than  m  ~  3.  A  study  of  these  high 
frequency  oscillations  will  be  reported  in  future 
papers. 


V.  Summary 

Preliminary  measurements  are  described  of 
electron  number  density  obtained  in  tbe  near 
field  of  a  200W  Busek  Hall.  The  approach  taken 
is  ultra-high  high  frequency  (3.3  mm 
wavelength,  90  GHz)  microwave  interferometry 
that  affords  high  spatial  resolution  suitable  for 
studying  the  near  exit  region,  where  the  plasma 
density  can  exceed  10llcm'3.  A  90  GHz  beam 
can  be  focused  down  to  a  waist  of  approximately 
7  mm,  and  over  a  50  mm  plasma  path  length,  a 
10Mcm'3  plasma  density  gives  rise  to  an  easily 
measurable  3°  phase  shift.  The  system  is  of  a 
phase-bridge  design,  utilizing  two  signal  arms 
split  from  a  fixed  frequency  source  (one  passing 
through  the  plasma)  that  recombine  at  two 
balanced  mixers.  A  line-integrated  electron 
density  is  obtained  by  comparing  the  two  signals 
from  the  mixers.  This  interferometer  is  shown 
to  be  suitable  for  measuring  time-dependent 
plasma  density  fluctuations  offering 
unprecedented  information  about  plasma 
fluctuations  in  the  near  exit  region,  and  the 
opportunity  to  study  turbulence-enhanced 
electron  transport  in  regions  where  the  plasma  is 
presumably  collisionless  and  free  of  interactions 
with  the  channel  wall. 
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Figure  10.  Fourier  decomposition  of  the  discharge  current  (top)  and 
microwave  mixer  signal  (bottom)  with  the  microwave  beam  passing 
through  the  near  field  (z  =  45  mm). 
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